home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 007 / threed.arc / COSX.F77 next >
Encoding:
Text File  |  1986-03-28  |  1.1 KB  |  52 lines
  1.     PROGRAM COSX
  2. C
  3. C Program to plot a three dimensional COS(X)/X function with
  4. C hidden-lines removed
  5. C
  6.     COMMON ANGA, ANGB, HV, D
  7.     COMMON DL, DM, DN, CX, CY, CZ, QX, QY, QZ, DD
  8. C
  9.     REAL X1, Y1, DX, DY, DELMIN
  10.     DATA X1, Y1, DX, DY /-16.0, -16.0, 0.5, 0.5/
  11.     INTEGER NX, NY
  12.     DATA NX, NY /65, 65/
  13.     REAL ADEG, BDEG
  14.     DATA ADEG, BDEG /60.0, 30.0/
  15. C
  16.     REAL A(100, 100)
  17. C
  18.     DELMIN = DX
  19.     IF(DELMIN.GT.DY) DELMIN = DY
  20.     WRITE(6,50)
  21. 50    FORMAT(1X,'Generating Cos(X)/X function values',/)
  22.     X = X1
  23.     DO 129 I = 1, NX
  24.         Y = Y1
  25.         DO 119 J = 1, NY
  26.         XSQR = SQRT(X * X + Y * Y)
  27.         IF(XSQR.LT.0.001) XSQR = 0.9 * DELMIN 
  28.         A(I,J) = COS(XSQR)/XSQR
  29.         Y = Y + DY
  30. 119        CONTINUE
  31.         X = X + DX
  32. 129    CONTINUE
  33. C
  34. 1    WRITE(6,100)
  35. 100    FORMAT(1X,'INPUT TILT AND ROTATION ANGLES ')
  36.     READ(5,*)BDEG,ADEG
  37. C200    FORMAT(2F10.2)
  38.     IF(BDEG.LE.90.0) GOTO 118
  39.     WRITE(6,101)CHAR(27)
  40. 101    FORMAT(1X,A1,'mA')
  41.     STOP
  42. C
  43. 118    PI = 3.141593
  44.     ANGA = ADEG * PI / 180.0
  45.     ANGB = BDEG * PI / 180.0
  46.     HV = 5.0
  47.     CALL THREED(A, NY, NY, 3)
  48.     CALL    PICSIZ(0.,0.)
  49.     GOTO 1
  50. C
  51.     END
  52.